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Abstract 



The formalism of jet modification in the higher twist approach is modified to describe a hard parton propagating 
through a hot thermalized medium. The leading order contribution to the transverse momentum broadening of a 
high energy (near on-shell) quark in a thermal medium is calculated. This involves a factorization of the perturbative 
process of scattering of the quark from the non-perturbative transport coefficient. An operator product expansion 
of the non-perturbative operator product which represents q is carried out and related via dispersion relations to the 
expectation of local operators. These local operators are then evaluated in quenched SU(2) lattice gauge theory. 

Keywords: Quark Gluon Plasma, Jet Quenching, Lattice Gauge Theory 



1. Introduction 

Jet quenching, the modification of hard jets fT) in a hot dense extended medium, is considered as one of the most 
penetrating probes of the Quark Gluon Plasma (QGP) formed in high energy heavy-ion collisions. At this time, the 
most promising approach to jet modification has been to treat the propagation of the hard patrons in the jet using 
perturbative QCD (pQCD), factorized from the soft medium which is treated non-perturbatively [2|. In so doing, the 
effect of the medium has been codified into a handful of transport coefficients: The transverse momentum squared 
gained per unit length, by a single parton, without radiation (q = {p])/t), the longitudinal momentum lost per unit 
length without radiation (e = {pL)lt), the fluctuation in longitudinal momentum per unit length (e2 - {p])lt) ||3]|4]|5l. 

As derived in the references mentioned above, all of these transport coefficients have well defined operator ex- 
pressions in a given gauge and also have gauge invariant definitions derived in effective field theory |6, 7|. However, 
to this day, there is no first principles calculation of these coefficients within QCD at the temperatures of relevance at 
RHIC and LHC collisions. In these proceedings, we report on such an attempt |8l. 

2. A single parton in a QGP brick 

In this attempt, the simplest example of jet modification in a finite thermaUzed medium will be considered: A 
single quark propagating through a box of length L, held at a fixed temperature, undergoing a single scattering which 
imparts transverse momentum to the quark. The quark then exits the medium and is observed. The length of the 
medium will be considered to be short enough such that the quark does not radiate during traversal. The coupling of 
the jet with the medium will be considered to be small enough such that secondary scattering is minimal. 

Consider a quark in a well defined momentum state |^)(= \q^ ,q^ , Oj_)) impinging on a medium \n) and then exiting 
with transverse momentum k_i_. 




(1) 
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The medium state absorbs this change in momentum and becomes \X}. The quark is assumed to be space-like ofF-shell 
with virtuality - 2q^q^ < with the negative z-axis defined as the direction of the propagating quark. The spin 
color averaged transition probability (or matrix element) for this process, in the interaction picture, is given as 

W(k) = ^ ^|!le-/'^"<<^-;„|r*e'/o'"^'«|^- + k^,X}{q- + k^,X\Te--£'"'''%-,n}, (2) 
^ /TV,. Z 

where, we have averaged over the initial color and spin of the quark. Assuming that the medium is in a thermalized 
state we have also averaged over the various initial states of the medium, weighted with a Boltzmann factor (Z 
represents the partition function of the medium without the incoming quark). Summing over the various values of k±, 
weighted with kj_, and dividing by the length yields the transverse momentum diffusion coefficient. 




Expanding the time evolution operator to leading order in coupling, converting factors of into transverse gra- 
dients (and finally into field strength operators), and summing over the final state \X}, we obtain the non-perturbative 
definition of the transport coefficient q, i.e., 

g'=-^j (2;r)3 -"^ {n\^^F ' ^{y )F ^(Q)\n). (4) 

In the above expression, there is no ordering between the two field strength operators. The expression above is 
not gauge invariant, but is gauge covariant. As a result, if one were to carry out an operator product expansion in 
terms of local operators, one could reorganize the expansion to only contain gauge invariant local operators. Any 
gauge dependence would then only be contained in the coefficient functions. The operators in the product above are, 
however, separated by a light like distance. As a result, a straightforward expansion in y^ will not suffice. 



3. Analytic Continuation and Finite Temperature Dispersion Relations 

The expression for q above may be written as the discontinuity of a more generalized operator product. 



q — Disc 



Disc 



N, 



r d yd . 



^J2q- 



(q + k)^ + ie 



(5) 



The generalized quantity Q should be considered as an analytic function of q^,q^. Keeping q^ constant as the large 
scale Q in the problem, we consider 2 as a function of q^. The jet transport coefficient is the discontinuity in Q in the 
region where \q^\ <K \q^\. 

We now consider the integral. 



dq^ Q(q^) 
Iniiq^ + QoY 



(6) 



where Qo is positive and of the order for . The contour is taken as a small counter-clockwise circle around the point 
q'^ = -Qq. The residue of this integral is given as. 



(7) 



When q^ = Qo ~ Q, the hard scale in the problem, one can expand the denominator in Q as a series of local operators 
involving ever higher powers of derivatives. 



NcQ^ 



-q ■ iD - Dl 



(8) 



/Nuclear Physics A 00 (2012) i-g] 



3 



We can now deform the contour to have it encircle the cut in the region where |^^| «c 2 the hard scale. The disconti- 
nuity in this region, is related to the discontinuity of Q which in turn leads to the physical transport coefficient q. In 
this region, however, one cannot use the expansion of Eq. (|8]l. Instead one obtains, 

^1 " ^1 Tr^^ ^~ / J ~^ , + ^ ^ , = dg + dg V{q ). (9) 

Nc J (Inr 2g „ ^ + 6o) J-a^-q g^ + Go Jo 

The second integral in the equation above, represents the vacuum process of a time-like quark decaying by the emis- 
sion of a hard gluon. The dimensionless quantity A in the limits of integration above is meant to indicate a small 
number A <s: 1. It merely states that the physical g integral is over a small range in g^, far from g^ = Qo- For a jet with 
maximum virtuality fi^ and (-) momentum g^, A^Q - ji^ jilg^). Given this small range of values of g^, the integrand 
can now be expanded as a Taylor expansion in g^ IQq. This series can then be compared with the series of Eq. ([HJ, 
and terms equated to obtain the various moments of g. 

The methodology outlined above can be made even more precise and straightforward by setting a definite value 
for 2° - g^ ■ While this will readjust the relative importance of the various terms in the series it allows for simpler set 
of operators that need to be evaluated numerically. This simplifies I\ in Eq. (|8]l to, 

^ n 11=0 ^ ' ' 

and similarly simplifies Eq. ([9]) with Qq replaced by g^ . For a virtuality //^ such that A^^^^ <k //^ <k )2, we can 
define a g^ or virtuality averaged g as. 



J-o* 



ag g\q ) - z.gi^ + 



g(Q")2Q^= dg^gig^) ^ 2gQ^ + ^ T , dD 



where the last partial equality is only valid in the limit that § is a slow function of g^, note = n^/g . One can now 
compare coefficients of g^ between Eq. ( [TO] i and Eq. (|9|, using the mean values defined in Eq. ([TTJ. 



4. Lattice calculation and Results 

The remaining methodology consists of simply evaluating the series of operator products in Eq. ( fTO] ) on an lattice 
at finite temperature and then comparing them with terms with the same power of g^ in Eq. ([9]). In order to calculate 
the thermal expectation values of the various local operators we have to simply rotate the time coordinate to the 
imaginary time direction: jc'^ — > ix*. A" — > /A"* and F^' iF'^'. 

For this exploratory study we calculate in quenched S U(2) gauge theory, as a simpler substitute for QCD. Calcu- 
lations are carried out on a lattice of spatial size = x 4, with temporal extent nj varying from 3 to 6. The scale 
is set using the renormalization group formula ll9l fT0l . 

In the equation above, ai is the lattice spacing, g represents the bare lattice coupling and A/, represents the one 
dimension-full parameter on the lattice. Comparing with the vacuum string tension, we have used A/, = 5.3 MeV. For 
a lattice at finite temperature or one with n, <K tis, the temperature is obtained as 

r=— . (13) 

n,aL 

The results reported in this proceedings are the result of 5000 heat bath sweeps per point. In the left panel of Fig.fT] 
we present the results of the lattice calculation for the first operator product in Eq. ( fTO] ), i.e., {F'^'F'^') scaled by 
to make it dimensionless. With these statistics, we find acceptable scaling with lattice size, especially in the region 
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above and away from the phase transition, i.e., for T > 400 MeV. In the right panel we calculate the second term in 
the expansion, {F^'TPF^'), in Eq. ( fTO] ), both with (green diamonds) and without (red diamonds) the large factor of 
in the denominator. While this operator product is by itself larger than the leading term, the large factor of in the 
denominator makes it considerably smaller than the leading term in the region with T < 600 MeV. As a result, in the 
region with 400 MeV< T < 600 MeV, the leading order term in Eq. ( [TO] i will yield an usable estimate for q. 
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Figure 1: Left panel: The temperature dependence of the local operator scaled by to make it dimensionless. Right panel: temperature 

dependence of absolute values of the the local operator (F^'F^') and the next-to-leading operator products ([F^'(2)''F+'j> [red diamond (dashed)] 
and (iF^'iiy^F^' Iq^) [green diamond (dot dashed)]. 



To estimate q we use the point at T = 400 MeV, where {F^'F^') = 0.01 GeV*. We are considering a lattice with 
a length given by 4 x n,ai - 4/0.4 GeV"' = 10 GeV"'. This states that the maximum virtuality of a jet (with a = 
20 GeV) which traverses such a length without undergoing radiation is given as //2 = EjL = 20/10/ V2 ^ 1.4 GeV^. 

Thus g+ = 1.4/40 GeV. We can now use the formula q = ^^^|^<F+'F+'). Using Q',(1.4GeV^) = 0.375 ifTTl . 
we obtain q - 0.134GeV^/fm for an SU{2) quark traversing a quenched SU{2) plasma. In most phenomenological 
estimates one quotes the q of the gluon. If the above calculation were done for an S U{2) gluon, the q would differ only 
by the overall Casimir factor of Ca/Cf = 2N^/(N^ - 1) = 8/3 yielding a §g = 0.358 GeV^/fm, at a T ~ 400 MeV 
(Note, this value is different from that in Ref. ID by a factor of 4, due to the neglect of this factor in that reference). 
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